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In the standard cosmological model a mysterious cold dark matter (CDM) component dominates 
the formation of structures. Numerical studies of the formation of CDM halos have produced 
several robust results that allow unique tests of the hierarchical clustering paradigm. Universal 
properties of halos, including their mass profiles and substructure properties are roughly consistent 
with observational data from the scales of dwarf galaxies to galaxy clusters. Resolving the fine 
grained structure of halos has enabled us to make predictions for ongoing and planned direct and 
indirect dark matter detection experiments. 

While simulations of pure CDM halos are now very accurate and in good agreement (recently 
claimed discrepancies are addressed in detail in this review), we are still unable to make robust, 
quantitative predictions about galaxy formation and about how the dark matter distribution changes 
in the process. Whilst discrepancies between observations and simulations have been the subject of 
much debate in the literature, galaxy formation and evolution needs to be understood in more detail 
in order to fully test the CDM paradigm. Whatever the true nature of the dark matter particle is, 
its clustering properties must not be too different from a cold neutralino like particle to maintain 
all the successes of the model in matching large scale structure data and the global properties of 
halos which are mostly in good agreement with observations. 
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1. INTRODUCTION: FROM COLD COLLAPSE 
TO HIERARCHICAL CLUSTERING 

1.1. A short history 

N-body simulations of the gravitational collapse of 
a collisionless system of particles pre-dates the CDM 
model. Early simulations in the 1960's studied the for- 
mation of elliptical galaxies from the collapse of a cold 
top-hat perturbation of stars 1 * 2 -! 3 -. The resulting virial- 
isation process gave rise to equilibrium structures with 
de Vaucouleurs 4 - or Einasto££ type density profiles, simi- 
lar to observations of elliptical galaxies. It is remarkable 
that the end state of almost any gravitational collapse, 
independent of the small scale structure and hierarchical 
merging pattern, leads to a similar global structure of the 
final equilibrium systemic. 

Computer simulations in the 70's attempted to follow 
the expansion and a collapse of a spherical overdensity to 
relate to the observed properties of virialised structures 
such as galaxy clusters 1 ^. Using a random distribution 
of particles with a Poisson power spectrum lead to the 
initial formation of many bound clumps, however it was 
observed that these bound structures were destroyed as 
the final system formed - resulting in a smooth distri- 
bution of matter. This overmerging problem persisted 
for over two decades and motivated the development of 
semi-analytical models for galaxy formation^ 1 -. 

During the 1980's, it was proposed that cosmic struc- 
ture formation follows a dominant, non-baryonic cold 
dark matter (CDM) component^ 2 -. Cold dark mat- 
ter could consist of new and yet undiscovered weakly- 
interacting massive particles (WIMPs), which occur for 
example in super-symmetric extensions of the Standard 
Model of particle physics 13 . "Cold" means that these 
particles have rather small thermal velocities, which al- 
lows the formation of very small structures, typically 
down to far below one solar mas o 14 i 15 ' 16 . CDM together 
with the even more mysterious dark energy (usually de- 
noted " A" ) are the dominant components of the ACDM 
model, in which all the ordinary matter accounts for only 
4.6 percent of the total. ACDM has by now become the 
"standard cosmological model" and its parameters (and 
therefore the initial conditions for structure formation) 
are now known to a reasonable precision^!. 

Computer simulations allow to follow the non-linear 
evolution of perturbations, starting from realistic and 
well constrained cosmological initial conditions. The fi- 
nal quasi-equilibrium structures are the dark matter ha- 
los that are observed to surround galaxies and galaxy 
clusters. During the 1980's, the first simulations of the 
CDM model were carried out. Large cubes of the universe 
were simulated in an attempt to match the large scale 
clustering of galaxies. Some of the most basic properties 
of collapsed structures were discovered - the distribution 
of halo shapes, spin parameters et o 18 ' 19 . It was not un- 
til the simulations of Dubinski & Carlberg that individ- 
ual objects were simulated at sufficiently high resolution 
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FIG. 1: Density profile of the million particle dark matter halo 
simulation of Dubinski & Carlberg 1990 (crosses). The solid 
line shows the best fit NFW profile (Eqn. [1} to the original 
data. This Figure was adapted frorr>22 by John Dubinski and 
it is reproduced here with his permission. 



to resolve their inner structure on scales that could be 
compared with observations 2 ^. Using a million particle 
simulation of a cluster mass halo run on a single work- 
station for an entire year, these authors found central 
cusps and density profiles with a continuously varying 
slope as a function of radius. They fit Hernquist profiles 
to their initial simulations but an NFW profile 21 provides 
an equally good fit (see Figure 1). Most likely due to a 
large softening length, the final virialised structure was 
almost completely smooth. 



Navarro et al. (1996) published results of simulations 
of halo density profiles from scales of galaxies to galaxy 
clusters. They demonstrated that all halos could be rea- 
sonably well fit by a simple function (Eqn. [I} with a con- 
centration parameter that was related to the halo mass 21 -. 
However, with less than 10 4 particles only the mass pro- 
file beyond about 5-10 percent of the virial radius was re- 
solved reliably. Shortly afterwards, simulations with 10 6 
particles showed cusps steeper than r~ l down to their 
innermost resolved point near one percent of the virial 
radius 23 . These simulations also resolved the overmerg- 
ing problem^ 4 . - the resolution was sufficient to resolve 
cusps in the progenitor halos enabling the structures to 
survive the merging hierarch y 23 i 25 i 26 . The final surviving 
substructure population is a relic of the entire merger his- 
tory of a given CDM halo. 
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1.2. State of the art simulations - convergence or 
discrepancies? 

Algorithmic and hardware development have increased 
the mass and spatial resolution by orders of magnitude 
(parallel computing, special purpose hardware, graphics 
pipelines etc). The first simulations used just a few hun- 
dred particles with length resolution that was a large 
fraction of the final structure. Today we can simulate in- 
dividual collapsed structures, in a full cosmological con- 
text with up to 10 9 particles and spatial resolution that 
is better than 0.1% of the virialised region and ~ 10 5 
substructure halos can be resolve d 27 ' 28 ! 29 . 

The first billion particle halo simulation "Via Lactea 
II" (VL-II hereafter)* 7 - was published 2008. About half a 
year later Springel et al. posted a preprint (0809.0898vl) 
in which significant discrepancies between VL-II and 
their Aquarius simulations were claimed. For some rea- 
sons these authors failed to acknowledge that the Aquar- 
ius and VL-II simulations are perfectly consistent, which 
has caused some confusion. Below we clarify the situation 
and to we hope to restore confidence in the robustness 
of the current state-of-the-art simulations like VL-II 27 , 
GHALO 28 and Aquarius* 9 -. 

Springel et al. 2008 claim two main discrepancies be- 
tween VL-II and Aquarius: subhalos in Aquarius are 
more concentrated than in VL-II and the abundance 
of subhalos as a function of their peak circular velocity 
(see Section f2 . 3 . 1 1) is about 30 percent higher in Aquar- 
ius. However, the Aquarius simulations still adopt the 
WMAP-lyr cosmological parameters, especially as = 
0.9, n s = 1.0 (see e.gi 17 ' 30 for definitions and current con- 
straints). Both values are significantly higher (about 3 
standard deviations) than the latest WMAP-5yr values* 7 - 
and introduce much higher typical amplitudes for the 
small scale fluctuations in the Aquarius initial conditions. 
VL-II on the other hand used the newer WMAP-3yr val- 
ues as = 0.74, n s = 0.95, which is somewhat lower, but 
consistent (within one a) with the latest parameters. The 
earlier VL-I simulation** had <jg = 0.74, n s ~ 0.9, due to 
an error in the initial conditions generator GRAFIC2 32 , 
which leads to somewhat lower (sub-)halo concentrations 
and velocity functions. 

The substantially different amplitudes of the typical 
initial small scale fluctuations have to be taken into ac- 
count when the properties of VL-II and Aquarius subha- 
los are compared. It is well known that halos and subha- 
los form earlier in the WMAP-lyr cosmology which leads 
to higher halo and subhalo concentrations**. The concen- 
tration difference between WMAP-lyr and 3-yr matches 
exactly with the "discrepancy" claimed by Springel et al. 
2008, i.e. the concentrations in VL-II and Aquarius dif- 
fer by precisely the expected amount. In other words, 
there is no discrepancy between the simulations, just 
a difference in adopted cosmological parameters. The 
older cosmology used in Aquarius simply leads to signifi- 
cantly higher halo and subhalo concentrations compared 
to what's expected in a ACDM Universe with up-to-date 



parameters**. 

Higher subhalo concentrations also cause a higher sub- 
halo peak circular velocity at the same subhalo mass, i.e. 
a higher subhalo velocity function is expected in a cos- 
mology with too much small scale power—. Therefore 
the higher subhalo velocity function, the second claimed 
'discrepancy', is also an expected consequence of the ex- 
cess small scale fluctuations in the cosmology chosen for 
the Aquarius simulations. 

Even if both simulations would have adopted the same 
cosmology, a 30 percent difference in subhalo velocity 
functions of individual host halos would not be a sig- 
nificant discrepancy, since velocity functions have sub- 
stantial halo-to-halo scatter—. The small scatter found 
among the six pre-selected Aquarius halos cannot be gen- 
eralised to the full population of galaxy halos, which 
shows much larger variation 3 ^. 

After allowing for the differences caused by the dif- 
ferent cosmologies the main results of Via Lactea and 
Aquarius, which used different initial condition gener- 
ators, different simulation codes and different analysis 
tools, are consistent. This is a reassuring confirmation 
of the robustness and maturity of cosmological N-body 
simulations. Current simulations arc excellent approxi- 
mations to cold dark matter halos on a large mass and 
spatial range. For most applications the uncertainties in 
N-body results (e.g. on the very inner density profiles, 
see Section 12.2.2)) are already much smaller than the ex- 
pected, but still poorly constrained effects of galaxy for- 
mation on the dark matter distribution, see Section [3] 



2. RESULTS FROM COLLISIONLESS 
SIMULATIONS 

Before we summarize the main simulation results on 
the properties of dark matter halos we introduce defini- 
tions of the mass and radius of dark matter halos. Decid- 
ing which material belongs to a halo and what lies beyond 
it is a non-trivial question, which has not received much 
attention until quite recentl y 36 i 37 i 38 i 39 . 

The usual halo definitions are still based on the sim- 
ple spherical collapse pictur e 40 ' 41 . These models assume 
that there is no kinetic energy in an overdense sphere 
of matter as it turns around and that this sphere virial- 
izes, rather than falling back radially into one point. The 
virial theorem predicts that the final halo radius is 0.5 of 
its turnaround radius and that this "virial radius" would 
enclose an overdensity of 178 times the critical density 
of the universe p cr it in an Einstein-deSitter cosmology 
(Qm = 1-0). For the now favored ACDM cosmology 
this overdensity is close to 100 times p cr it or roughly 300 
times the mean matter density pm- However CDM ha- 
los form very differently than assumed in the classical 
spherical collapse model: material ending up near this so 
called "virial radius" did not undergo a collapse by any- 
where near a factor of two 3 - 7 -. The ratio of final radius 
to turnaround radius is much larger in the outer halo 



FIG. 2: Dark matter density maps from the "Via Lactea II" (VL-II) simulation 27 . Cubes of 800 proper kpc are shown at 
different redshifts, always centered on the main progenitor halo. VL-II has a mass resolution of 4,100 Mq and a force resolution 
of 40 pc. Initial conditions were generated with a modified, parallel version of GRAFIC322. The high resolution region (some 
of its border is visible in the upper panels) is embedded within a large periodic box (40 comoving Mpc) to account for the large 
scale tidal forces. More images, movies and data are available at http://www.ucolick.org/~diemand/vl/ 



and typical orbits extend out to 90% of the turnaround 
radius—, i.e. well beyond the formal virial radius. All the 
conventional overdensity definitions, by which the halo 
radius encloses of the order of 200 times p C rit or pu, are 
therefore too small, and they underestimate the extent 
and mass of dark matter halos significantly^. However, 



spherical overdensities are easy to measure in simulations 
and useful for comparisons. In the following we adopt the 
largest of these widely used halo radius definitions: r2oo, 
defined to enclose 200 times pm- As a proxy for halo 
mass we will use M 2 qo — M(< r 2 oo)- 

One has to keep in mind that practical, ad hoc halo 
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FIG. 3: Dark matter density map within the inner 200 kpc of 
the "GHALO" simulation^. This galaxy scale halo is resolved 
with over one billion particles of 1000 Mq each. It contains 
over 100,000 subhalos, the largest are visible as bright spots 
in this image. 

definitions like r2oo and M200 are rather different from 
the concept of halo mass defined in theoretical mod- 
els of structure formation. This can be confusing and 
has affected many attempts for comparisons of theory 
with simulations, e.g.: halo mass functions (see 12. ip : 
halo mergers and merger trees are well defined in the- 
ory, but not in practice since the M200 of the remnant 
may be smaller than the sum of its progenitors M 2 oo 
and because some (sub-)halos do not stay within r2oa, 
i.e. they de-merg o 37 ' 39 ' 42 ' 43 ; by convention we will also 
refer only to halos within the virial radius of a larger 
halo as " subhalos" , but this actually miss-classifies many 
true satellite halos beyond this ad hoc scale as field halos, 
which introduces environmental dependencies into the to- 
tal "field" halo population 37 ; the halo mass M2oq(z) and 
radius r2oo (z) grow non-stop by definition, even when no 
mass is accreted^. 



2.1. Mass function of halos 

For a range of practical halo mass definitions, e.g. 
based on a given spherical overdensity, the abundance 
of field halos as a function of mass and redshift can be 
accurately measured in cosmological simulations of large, 
representative, periodic volume a 44 ' 45 . 

The classic analytical approach by Press & Schechter— 
combines the statistics of a hierarchical, Gaussian ran- 
dom field with the spherical, radial top-hat collapse 
model and it predicts the abundance of collapsed ob- 
jects as a function of mass. This idealized model matches 



the halo abundances measured in simulations surprisingly 
well, but it predicts too many small halos (M < M*, 
where M* is the mass of a typical, one-sigma halo form- 
ing now). This difference is usually interpreted as being 
caused by the assumption of a spherical collapse. Allow- 
ing for an ellipsoidal collapse introduces free parameters, 
which allows for a very accurate fit to the measured mass 
function o 44 ' 47 . 

Alternatively, one can argue that the main difference 
between the model and the simulation is the collapse (or 
lack thereof): the infalling mass hardly loses any energy, 
but typically orbits back out to 90% of its turnaround 
radius 39 , i.e. significantly beyond the "virial" radius. In 
other words, one cannot expect to find the collapsed mass 
predicted by the Press-Schechter approach inside the 
overdensity radii commonly used to define halo masses 
in simulations. A recently suggested improvement is to 
use the stationary mass, i.e. the mass within the largest 
radius with zero mean radial velocit y 36 ' 38 . The static 
mass exceeds the "virial" mass substantially in halos be- 
low M* and at z=0 the Press-Schechter mass function 
fits the abundance of halos as a function of their static 
mass very well 3 - 8 -. At z=l and z=2 the agreement is less 
impressive, the Press-Schechter mass function lies above 
simulated abundances 3 - 8 -. 

Further improvements in these directions seem possi- 
ble, since Press-Schechter predicts the collapsed mass, 
not the static mass. Above M„ accretion typically sets 
in near r v i r , i.e. static and virial radius are similar 36 . But 
also above M* significant amounts of previously collapsed 
material are found beyond the virial radiu o 37 ' 39 ' 42 ' 43 , so 
both the static and the virial mass seem to underestimate 
the collapsed mass for halos above M*. 



2.2. Global halo properties 

Here we briefly summarize some basic, global proper- 
ties of CDM halos: concentration and density profiles; 
shapes, spin and velocity distribution. 

2.2.1. Halo formation, density profiles and concentrations 

Self-similar infall models predict scale free, nearly 
isothermal (p cx r -2 ) profile d 40 ' 41 ' 48 ' 49 . Simulated pro- 
files are found to be steeper than p oc r~ 2 in the outer 
parts and shallower in the inner regions 20 , i.e. their cir- 
cular velocity profiles V c (r) = ^GM{< r)jr have a well 
defined peak, which serves as a natural halo size scale 
Vmax and scale radius rymax- 

Self similar spherical infall maintains a distinct rela- 
tion between the initial conditions and the final structure, 
while violent relaxation assumes information on the ini- 
tial conditions is chaotically lost. Simulations show that 
despite their non-isothermal density profiles, CDM ha- 
los still do have much in common with idealized spheri- 
cal infall models: i) particles that collapse in rare, early 
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peaks (e.g. formation sites of first stars and old globu- 
lar clusters) are located closer to the potential minima 
of the entire turnaround region due to peaks biasing and 
those particles do end up closer to the center of the final 
halo 50 ' 51 , ii) the typical particle apocenter distances are 
close to their turnaround radii 3 -^. Modified infall models 
are indeed able to reproduce some of the features of halo 
density profiles found in cosmological simulations 5 ?. 

Over a mass range spanning 20 decades, from micro- 
halos to galaxy clusters, the spherically averaged CDM 
halo density profile can be approximated with the same 
universal form (NFW—): 



p(r) 



(r/r s )T(l + r/r s ) 3 "T : 



(1) 



where 7=1. The scale radius r s is related to the peak 
circular velocity scale by fvmax — 2.163 r s and it is used 
to define the halo concentration c v i r = r v ; r /r s , where 
the "virial" radius r v ; r is defined following one of the 
ad hoc overdensity criteria described above. Halo con- 
centrations Cvir (and equivalently scale densities p s ) are 
related to the halo formation time: early forming halos 
tend to have higher c v i r and p s at z=0. In CDM the typ- 
ical amplitudes of density fluctuations a(M) decreases 
from dwarf to galaxy cluster scales, i.e. smaller halos 
form earlier on average than larger ones and they end 
up having higher median concentration a 21 ' 33 ' 53 ' 54 ' 55 . At 
a given mass, the concentrations of individual halos have 
a large scatter: the standard deviation in log c v i r is 0.18. 
On subsolar mass scales the CDM power spectrum ap- 
proaches P(k) cx k~ 3 , i.e. a(M) approaches a constant, 
which leads to very similar halo formation times and 
halo concentrations over a wide range of masse s 14 ' 56 ' 57 . 
Even Earth mass micro-halos, the first and smallest sys- 
tems in the CDM hierarchy, have NFW-like density pro- 
files. A systematic study of their typical concentrations 
is still lacking, but values found in the small sample of 56 
(c v i r (z — 0) ~ 80) seem consistent with the predictions 
of the Bullock et al. model 57 . 

A simpler and more general measure of halo concentra- 
tions is the mean density within rv mm . It is well defined 
both for isolated halos and subhalos and it is independent 
of assumptions on their "virial" radius or their density 
profile 3 - 7 -: 



cv 



_ P (< TVmax) 
Pcrit,0 



v m 



47rGpc.it, 



(2) 



For an NFW halo it is easy too convert from cy to c V i r 
f 3 - 7 -] . Since the NFW form is not a very good fit to most 
CDM halos, the measured c v i r depend somewhat on the 
details of the fitting procedur e 33 ' 53 ' 54 ' 55 . These compli- 
cations could be avoided by using cy, which is a robust 
concentration measure as long as V max and fvmax & r e 
resolved in the simulation. 



2.2.2. Inner density profiles 

The NFW form (Eqn. [I} was proposed as a fit to 
CDM density profiles in the radial range from 0.01 r2oo 
to r2oo [ 21 ] ■ As larger simulations started to resolve scales 
around 0.01 r2oo it became clear that most halos are 
significantly denser in their inner parts than their best 
fit NFW pro file2 3 .i25 i S8i 5 ^^iS 1 . i 62 i 63 i 6i 1 65. Samples of high 
resolution, relaxed, isolated CDM halos also revealed sig- 
nificant halo-to-halo scatte r 64 ' 65 : some halos follow the 
NFW form quite well, while others are better approx- 
imated with a steeper profile as suggested by-2 3 -. Most 
halos lie somewhere in between and their average devia- 
tions from both fitting functions are larger than 10 per- 
cent (more than 20 percent at some radii)^ 5 -. Improved 
fits require an additional free parameter to account for 
the substantial halo-to-halo scatter. Good fits are ob- 
tained by letting the inner slope 7 of the NFW function 
(Eqn. [I} become a free parameter, instead of forcing 
7=1. Another option is to use the Einasto profile^i 6 -! 6 ^ 



p(r) = p s e 



W(r/r 8 )«-l] 



(3) 



where a is the additional free parameter. Down to 0.03 
rvmax both forms are very similar and both fit simulated 
CDM density profiles very wel&SiiS 5 .^. The two func- 
tions only differ significantly at very small radii (below 
about 0.6 percent of rv max , i-e. below about 500 pc in a 
galaxy halo). Only very few (if any) simulations are cur- 
rently able to resolve such small scales. The halos i n 27 ' 67 
are denser than their Einasto fits in the inner parts and 
they are well approximated by 7 ~ 1.2 cusps. The same 
is true for the billion particle GHALO down to 400 po^ 8 .. 
However, the higher mass resolution of GHALO might al- 
low significantly smaller scales to be resolved and a new 
functional form for the fitting function of the density pro- 
file has been proposed 28 : 



p(r) = p e- x ^ (1+r / R ^ 2 



(4) 



This function has a constant logarithmic slope down to 
a scale R\, beyond which it approaches the central max- 
imum density po as r — ► 0. If one makes a plot of 
d\np/dln(l + r/R\) vs. ln(l + r/R\), then this pro- 
file forms an exact straight line with slope — 2A. Good 
fits are obtained with A ~ 0.10 28 . 

However, to demonstrate convergence on such ex- 
tremely small scales, several issues have to be addressed 
first: One necessary condition is to resolve the short 
dynamical times in these inner, very high density re- 
gions. The widely used, but unphysical, time step cri- 
terion At cx y/e/\a\, where e is the force softening, \a\ 
the local acceleration, can lead to artificially low central 
densities^ 7 . The VL-II and GHALO simulations use an 
improved criterion based on the local dynamical time 68 . 
Another caveat is that the material ending up in the in- 
ner 500 pc of a galaxy halo comes mostly from very early 
forming, dense progenitor halos (4cr-peaks and higher—). 
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At the starting redshifts of these simulations the varia- 
tion on the smallest resolved mass scales a(M m [ n ), is in 
the range 0.1 to 0.2. 4er-peaks already have the (mildly) 
nonlinear density contrast of 0.4 to 0.8, which artifi- 
cially lowers their formation rcdshift and scale density. 
The common practice of starting the simulations when 
tr(M m i n ) ~ 0.1 artificially lowers the halo abundace at 
high redshifts^ 5 . Furthermore, relaxation effects will be 
present within these first structures to collapse since they 
are always resolved by a small number of particle s 69 ' 70 . 

Beyond 0.6 percent of rymax there is a general con- 
sensus that halos are denser than NFW and that they 
are well fitted by functions which are steeper and cus- 
pier than NFW at these radii like Eqn. [T]with 7 ~ 1.2, 
or the Einasto form (Eqn. [3]) with 7 ~ 0.17. To resolve 
the remaining uncertainties in the density profiles of pure 
CDM halos below 0.6 percent of rv max is of some numer- 
ical and theoretical interest, but it does not affect the 
observational predictions from CDM significantly: 

• In galaxies, groups and clusters these scales are 
dominated by the baryons. Predictions on the dark 
matter and total mass distribution require a real- 
istic treatment of the baryons and their dynamical 
interactions with the dark matter (see Section [3]). 

• The predicted kinematics of dark matter dominated 
dwarf galaxies are indistinguishable for these two 
density profiles, even if a large numbers of accurate 
proper motion measurements become available 7 ^. 

• The annihilation properties of small dark matter 
dominated halos (see Section 14. 1|) do not change 
significantly: The resulting annihilation luminosi- 
ties of both fits and of the measured profile in Fig- 
ure |4] all lie within 5 percent and the half light radii 
are within 10 percent of each other. The best NFW 
fit to the VL-II density profile however underesti- 
mates the halo luminosity by a factor of 1.4 and it 
overestimates the half light radius by 1.3. 




FIG. 4: Density profiles of Via Lactea halo and several sub- 
halos. Main panel: Profile of the main cold dark matter 
halo (thick line) and of eight large subhalos (thin lines). The 
lower panel gives the relative differences between the simu- 
lated main halo profile and the Einasto fitting formula (Eqn. 
[3]), with best fit parameters: a — 0.170, r 3 = 21.5 kpc, 
p s = 1.73 x 10 -3 M pc -3 (red curve) and function (flj 
with a best fit inner slope of 7 = 1.24, r s = 28.1 kpc, 
p s = 3.50 x lO" 3 M pc" 3 (blue curve). The vertical dot- 
ted line indicates the estimated convergence radius of 380 pc: 
simulated local densities are only lower limits inside of 380 
pc and they should be correct to within 10% outside this re- 
gion. The cuspy profile is a good fit to the inner halo, while 
the Einasto profile has a too shallow slope in the inner few 
kpc, causing it to overestimate densities around 4 kpc and 
to underestimate them at all radii smaller than 1 kpc. The 
same behavior in the inner few kpc is found at also at higher 
redshifts, while the large residuals in the outer halos on the 
other hand are transient features. Inset: Rescaled host (thick 
line) and subhalo (thin lines) density profiles multiplied by 
radius square to reduce the vertical range of the figure. This 
figure was reproduced from2£. 



2.2.3. Shapes, spin and velocity distributions 

Another main result from cosmological collision-less 
simulations is that CDM halos are rather elongated, sig- 
nificantly more than for example elliptical galaxies 20 . 
The mean minor-to-major axis ratio is well described by 
< s >= 0.54(M v i r /M*) -0 050 , where s is measured at 
0.3 r v j r . The rms-scatter around the mean is about 0.1 
p 2 -] . The length of the intermediate axis is usually closer 
to the minor axis, i.e. most halos are prolate. Towards 
the center CDM halos are even more elongated, i.e. the 
axis ratios become smaller while their orientations remain 
fairly well aligned. At a growing radius like 0.3r v i r , halo 
samples become rounder with tim o 72 ' 73 , while individual 
halos have quite stable shapes at some fixed, inner radius 
(except of course during major mergers) 7 -^. 



CDM halos are supported by nearly isotropic velocity 
dispersions, not by rotation. The tangential and radial 
velocity dispersions measured in spherical bins define the 
anisotropy parameter (3 = 1- 0.5a^ ngcnti Jcr^ dilil . The 
typical, relaxed CDM halo has (3(r) ~ 0.35(r/r vir ) - 35 , 
i.e. (3 changes from nearly isotropic near the center to 
mildly radial in the outer halo 75 ' 76 . 

The velocity distributions of CDM halos are Gaussian 
only near the scale radius, i.e. where the density profile 
is roughly isothermal (p cx r~ 2 ). Analytical calculations 
predict more peaked distributions for the inner halo, and 
broader distributions beyond the scale radiu s 22 ' 76 ' 77 ' 78 
in agreement with the measured distributions in CDM 
simulations^* 7 ^. 

Halo shapes are supported by velocity dispersion el- 
lipsoids, which are elongated in the same direction as 
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the mass distribution, both globally and also locally 81 . 
Fairly symmetric velocity distributions in the tangen- 
tial directions are found throughout the hale 81 , i.e. 
there is about as much negative as positive angular mo- 
mentum material relative to any given reference axis. 
Small asymmetries lead to a some residual spin A = 
iJWE^/iGM 5 / 2 ), which has a median of only A = 0.04 
and a lognormal distribution of width o\ ~ 0.56, both 
for the dark matter and for adiabatic gas. However, the 
two spins are often poorly aligned, with a median mis- 
alignment angle of 30 degrees^ 2 -. The dark matter an- 
gular momentum tends to align roughly with the minor 
axis of the halo shape, with a mean misalignment of 25 
degrees 8 ^, while disk galaxy orientation and halo shape 
beyond 0.1 r V j r are completely uncorrelated 8 -^. It is of- 
ten assumed that the net halo angular momentum corre- 
lates with disk galaxy size and orientation, however only 
some selected fraction of the halo material with its wide 
variety of angular momenta can be incorporated into a 
realistic disk 82 . In our Galaxy for example, a signifi- 
cant fraction of the total available baryonic angular mo- 
mentum is not in the disk but in the polar orbit of the 
Magellanic clouds given their relatively large distances, 
masses and proper motions^. How exactly disk galaxies 
form out the angular momentum distributions available 
in ACDM halos remains an open question despite much 
recent progres o 86 ' 87 . 



2.3. Substructure 

A major contribution of N-body simulations to our un- 
derstanding of structure formation was to demonstrate 
how hierarchical merging gives rise to a vast amount of 
surviving substructure, both gravitationally bound (sub- 
alos) and unbound (streams). 

2.3.1. Subhalo abundance: velocity and mass functions 

The abundance of subhalos within r2oo = 402 kpc of 
in the VL-II simulation (Figure [5|) is well approximated 
by 

N{> K iax ) = 0.036 (V max /K iax , host )- 3 , (5) 

where Vinax.host = 201kms _1 . It is close to the median 
abundance found in a large sample of halos simulated 
using similar cosmological parameters^. In earlier halo 
samples using a higher normalization of the power spec- 
trum (as — 1.0) a higher normalization of 0.042 for the 
median abundance was found 8 - 8 -. This difference is con- 
sistent with the expected cosmology dependence 3 ^. Both 
samples demonstrate significant halo-to-halo scatter of 
about a factor of two. Some of the variation comes from 
the scatter in halo concentration, which introduces scat- 
ter in Ki iaXj host at a given host halo mass. Normaliz- 
ing the subhalo function to V r c (r 2 oo) (he. to M 2 oo) re- 
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FIG. 5: Subhalo and sub-subhalo abundances in the VL-II 
halo. Number of subhalos above Vinax within r2oo = 402 kpc 
(thick solid lines) and within 100 and 50 kpc of the galactic 
center (thin solid lines). The dotted line is N(> V max ) = 
0.036 (Fmax/Vmax.host) -3 , where Fmax.host = 201 kms -1 (at 
^Vmax.host = 60 kpc). It fits the subhalo abundance above 
F m ax — 3.5 km s -1 . The number of smaller subhalos is ar- 
tificially reduced by numerical limitations. The inset shows 
the sub-subhalo abundance within riooo (enclosing 1000 times 
the mean matter density) of the centers of eight (same ones 
as in Fig. [4} large subhalos (thin solid lines), riooo is well 
inside of the tidal radius for these systems. The thick solid 
line shows the subhalo abundace of the host halo inside of 
its riooo = 213 kpc. The (sub-)subhalo Vmax values are 
given in units of Viooo = \JGM(< riooo)/riooo of the cor- 
responding host (sub-)halo. Lines stop at Vmax = 2kms _ . 
The mean sub-substructure abundance is consistent with the 
scaled down version of main halo, and both the mean abun- 
dance and the scatter are similar to the values found i n 35 ' 88 
for distinct field halos. This figure is from2£. 

duces the scatter to about a factor of 1.25 [2sL[. With- 
out normalizing to the host halo size, the differential 
V„iax-distribution function of subhalos in dwarf galaxies 
to cluster halos scatters around 

dn/dV max = 1.5 x 10 8 T/" a 4 x 5 (/i 3 Mpc^ 3 kms" 1 ) , (6) 

again with a halo-to-halo scatter of about a factor of 
two2S. The median subhalo abundance given by the nor- 
malized ([5]) or non-normalized velocity function ^ seems 
to be approximately self-similar, i.e. independent of host 
halo mass and redshift 88 . Even tiny host halos (< Mq) 
show similar abundances, perhaps surprisingly so given 
that the density contrasts between subhalos and the host 
are much smaller on these scales due to their similar for- 
mation times caused by the nearly flat a(M ) 89 i 90 . 

While cosmological simulations are able to resolve the 
subhalo content of a given dark matter halo accurately, 
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the exact abundance of substructure around a given 
Galaxy, like for example the Milky Way, remains uncer- 
tain: 

• The Vmax.host of the CDM halo in which a 
galaxy with a rotation speed of 220 kms -1 would 
form could lie anywhere between 160-220 km s -1 
j9i 1 92 1 93] < The resulting unC ertainty in N(> V msix ) 

spans about a factor of (220/160) 3 = 2.6. 

• At a given Vmax.host, and also at a given M200, the 
subhalo abundance within r2oo has a substantial 
halo-to-halo scatte r 35 ' 88 . 

These theoretical uncertainties have to be considered 
when CDM subhalos are compared to the subhalos 
around the Milky Way satellite galaxies. 

Cumulative mass functions of subhalos can be ap- 
proximated by power-laws of the form M~ a , with a — 
1.9 to 2.0 and normalized so that the mass in subha- 
los larger than 10 _6 Mh os t is between about 5 to 15 

p Crcc m -23, 25. 26, 29, 31, 37,58, 75, 89, 94, 95, 96, 96, 97, 98, 99 rj^g s ^ ee p 

slope of the mass function means that there is a signifi- 
cant amount of mass in small subhalos, which are still un- 
resolved in current simulations. Numerical convergence 
studies show that about four hundred particles per sub- 
halo are required to resolve the mass function within T2oo 
f 7 - 5 -]. The convergence test in Springel et al. (2008) con- 
firm these results, however subhalos resolved with only 
60 particles are still included in their mass function*, 
which artificially lowers the slope a of their best fitting 
power-law. 

2.3.2. Subhalo evolution and their final spatial distribution 

The abundance of field halos of a given (moderate) 
mass is proportional to the dark matter density of the 
environment. This proportionality is altered as matter 
and halos fall into a host halo and tidal forces reduce the 
mass of systems, especially those orbiting close to the 
host halo center. Subhalos selected by their final remain- 
ing mass are therefore more extended than the matter 
distribution (see Figure ^SL^SUmm^iSLm. The 
number density of subhalos is independent of the mass 
threshold and roughly proportional to the dark density 
times radius: n sub , M0 (r) cx pdm(t) x r ^2£L2&22.] . V max 
is less affected by tide o 37 ' 103 and V max selected subhalos 
are show less of a difference relative to the dark matter 
distribution (see Figure [H]). 

During tidal mass loss rymax becomes smaller and 
the enclosed mean subhalo density cy increases. Sub- 
halos near the center of the host halo end up having 



* In the dense inner regions of halos and subhalos the numeri- 
cal requirement are actually even highesSSiS. In these regions 
however Springel et al™. include systems resolved with only 20 
particles in their sub-substructure abundance estimates. 



much larger concentrations cy than field halos (see Fig- 
ure 

In terms of a subhalos dark matter annihilation lumi- 
nosity L oc V^axAVmax °c ^maxV^v 7 (see Section [47TJ) , 
the increase in cy partially compensates for the reduced 
Vmax- The total annihilation luminosity from the entire 
CDM subhalo population therefore traces the dark mat- 
ter distribution, even tough the mass in subhalos is much 
more extended (see Figure |6|) . In other words the total 
subhalo luminosity mostly (except close to the halo cen- 
ter) follows the dark matter distribution, because L is not 
significantly affected by tidal stripping for most subhalos. 
This is not surprising, given that half of a subhalos anni- 
hilation originates from the rather small radius of about 
0.07 rymax- Smaller halos than those resolved in cur- 
rent simulations have higher typical concentrations, i.e. 
higher densities inside 0.07 rv max and their annihilation 
luminosity is expected to resist tidal losses even better. 

Tidal stripping removes mass from the outer, loosely 
bound regions of subhalo o 25 ' 31 ' 103 . Stripping is well ap- 
proximated by removing the mass beyond the tidal radius 
over some timescal o 34 ' 37 ' 104 . The tidal radius is defined 
so that the host halo density is equal (or similar to) the 
subhalo density at the subhalos tidal radius. High den- 
sity parts of subhalos are therefore able to survive in- 
tact even close to the center of the host. This explains 
why practically all (97% since z=l) subhalos survive until 
the present time, despite substantial mass loss in some 
cases2I. Subhalos have an inner, tightly bound region 
unaffected by mass loss, whose extent depends on the 
subhalo concentration and on its orbi t 31 ' 105 ' 106 ' 107 . The 
mass profile and the substructure content deep inside 
subhalos are the same as found in the inner parts of field 
halos^i. These two findings are related, because tidal 
stripping removes smooth and clumpy material without 
preference for one or the other, but from the outer parts 
only. The inner regions of subhalos retain their cuspy 
density profiles (Figure \^ 2 122£L i 105 i wr and their sub- 
structure (Figure [Spl. The impression that subhalos 
should have shallower inner density profiles 1 ^ and less 
substructure 29 than the inner parts of field halos is caused 
by insufficient numerical resolution. 

Subhalos move on rather radial orbits. The median 
peri- to apocenter ratio is 1:6 [ 25 ' 37 ], and only 5% of 
the orbits are rounder than 2:3 [21]. Only a few sub- 
halos are massive enough to suffer significant dynamical 
friction, which causes decaying orbits, disproportionally 
large mass loss and even complete merging with the cen- 
ter of the main halo in some cases2I. Most subhalos, 
and dark matter particles, move on fairly regular orbits: 
They reach nearly constant median apocentric distances, 
which lie close to their initial turnaround radii2 9 -. A few 
subhalos even gain energy during their pericenter passage 
in three-body interactions involving a larger subhalo and 
the host haloiS 9 .. 
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FIG. 6: Abundance and concentrations of subhalos vs. dis- 
tance from the galactic center in the VL-II simulation. Top: 
The number density profile of subhalos (symbols) is more ex- 
tended than the dark matter density profile p(r) (thick line). 
The difference relative to the matter distribution is largest 
for mass selected subhalo samples: Squares show a subhalo 
sample with a mass above 8 x 10 M©. Selecting subhalos 
by their peak circular velocity V ma x reduces the difference: 
Triangles are for subhalos with V max > Skms" 1 . Select- 
ing by subhalo annihilation luminosity (L oc V max Avmax > 
(5kms _1 ) 4 /(0.15kpc)) eliminates the difference throughout 
most of the halo (circles). Subhalo luminosity closely traces 
the dark matter distribution, except in the very inner regions 
(about 5% of r2oo, or 20 kpc for the VL-II halo). Bottom: 
Subhalo concentrations cv (median and 68% range are shown) 
increase towards the center, where the stronger tidal force re- 
move more of the outer, low density parts from the subhalos. 
To make sure their cy axe resolved, only subhalos larger than 
V m ax = 511ms" 1 are used. The error bars indicate the sta- 
tistical uncertainties in both panels. This figure was adapted 
from—. 



2.3.3. Subhalo shapes and orientations 

The shapes of subhalos are similar to those of field ha- 
los (see l2.2.5)l . but subhalos tend to be a bit rounder, es- 
pecially the ones near the host halo center—. Tidal inter- 
actions make individual subhalos rounder over timei^, 
and they also tend to align their major axis towards the 
center of the host halo 74 ' 111 ' 112 ' 113 . The alignment is of- 
ten maintained over most of the subhalos orbits, except 
during pericenter passages 74 . A similar radial alignment 
has been found for red galaxies in SDSS groups 114 . The 
major axes of the Milky Way dwarf satellites might also 
be preferentially aligned radially, i.e. roughly towards us. 




FIG. 7: Local phase space densities calculated with EnBiD 
in a 40 kpc cube in the center of the VL-II halo at z~0 (com- 
pare to the last panel in Figure [5}. Subhalos have very high 
central phase-space densities (> 10 -5 M pc -3 km -3 s 3 ) due 
to their steep inner density cusps and their relatively small 
internal velocity dispersions. The most clearly visible streams 
still have quite low densities (about 100 times below the lo- 
cal density) but owing to their low velocity dispersion (about 
10 times smaller than that of background particles) they just 
barely manage to stand out in local phase space density (these 
streams have about 10 -9 M© pc -3 km -3 s 3 ). 



If that is case current mass estimates based on spherical 
models (e.gi 115 i 116 i 117 i 118 ) would be biased towards some- 
what higher values. 



2.3.4. Other halo substructure: caustics, streams and voids 

Besides the gravitationally bound, dense subhalos dis- 
cussed above, there exists additional structure in the 
phase space of CDM halos. Current simulations are now 
starting to resolve some of this structure (see Figure [7]), 
although finite mass resolution and artificial numerical 
heatin g 69 i 70 still severely limits our ability to detect and 
resolve fine grained phase space structure. The coherent 
elongated features in Figure [7J are dark matter streams 
which form out of material removed from accreted and 
disrupted subhalos. In cases where the disrupted subhalo 
hosted a luminous satellite galaxy, the resulting streams 
would contain not only dark matter but also stars and 
produce detectable features in galactic stellar halo*^. 
This process explains the origin of stellar streams ob- 
served around the Milky Wa y 120 i 121 . 

Cosmological infall into a halo does at some times 
lead to the characteristic patterns in the f(v r , r) planed, 
which are predicted in the classic secondary infall 
model s 48 ' 49 . However, in real space significant density 
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enhancements (infall caustics) do not occur, since large 
random motions among the clumpy, infalling material 
prevent their formation 3 - 9 -. Cold infall is expected to oc- 
cur only in the first halos forming at the bottom of the 
CDM hierarchy The resulting caustics might propagate 
into larger halos at lower redshift*^. 



3. EFFECTS OF BARYONS 

There is a lot of work remaining to be done to quan- 
tify the effects that baryons can have in modifying the 
distribution of dark matter. Simulations become more 
complex and expensive and we do not as yet have a clear 
understanding of how galaxies form and evolve. 



3.1. Adiabatically steepening the density profile 

The dark matter density profiles can steepen 
through the adiabatic contraction due to dissipating 
baryon o 124 ' 125 ' 126 ! 127 . The strength of this effect depends 
on the baryonic fraction that slowly dissipates via radia- 
tive cooling. However, accretion of baryons via cold flows 
may dominate the growth of many galaxies 1 ^, thus it is 
not yet clear how strongly this changes the inner distri- 
bution of dark matter in galaxies. For a halo that cools 
the cosmologically available baryons into a disk compo- 
nent, the dark matter density at a few percent of the 
virial radius increases by about a factor of two and the 
final, total (baryons and dark matter) density profile can 
resemble an isothermal sphere - comparable to observed 
constraints on elliptical galaxies 1 ^. For isolated galax- 
ies less massive than the Milky Way, the baryon fraction 
decreases rapidly, Mb aryon oc V^ ir , such that the smallest 
galaxies have captured and cooled less than 10% of the 
available baryons 1 ^ 9 -. 



3.2. Dynamically flattening the central cusp 

At smaller radii, the growth of supermassive black 
holes or central nuclei can steepen or shallow the very 
central dark matter cusp depending of whether these 
structures grow adiabatically or through mergers. Gon- 
dolo and Silk (1999) 130 explored the effects of slow cen- 
tral black hole formation on the CDM cusp in the con- 
text of an enhanced annihilation signal. This mechanism 
can create isothermal cusps (p oc r~ 2 ) on parsec scales, 
with a boost factor of several orders of magnitude. If 
a dense star cluster surrounds the black hole (as in the 
center of the Milky Way) a shallower cusp (p oc r~ 1,5 ) 
is expecte d 131 ' 132 . On the other hand, if supermassive 
black holes form via mergers, then dark matter particles 
(e.g. WIMPs) can be ejected from the central halo region 
via three body encounters. This process could produce 
a constant density region of dark matter with a mass 



deficit roughly proportional to the mass of the sinking bi- 
nary objects^ 3 - 3 -. Similar behavior would result from the 
formation of central stellar nuclei in galaxies. Dissipa- 
tive growth would increase the central dark matter den- 
sity, but formation via merging of existing star clusters 
would lead to an inner corei 3 -^. A similar mechanism was 
studied in the context of cluster halos, whereby energy 
transfer to the dark matter background from dynami- 
cal friction acting on massive satellite galaxies, gave rise 
to an inner region with constant dark matter densityi 3 -^. 
All of these processes have yet to be studied in a realistic 
cosmological context. 



3.3. Feedback and stirring - flattening the cusp 

Feedback from the star formation process has fre- 
quently been invoked to flatten cusps, especially in 
dwarf galaxies which have challenged the CDM paradigm 
through observations of rotation curves, stellar velocities 
and star-cluster kinematics. A single violent event, which 
somewhat unrealistically ejects a cosmological baryon 
fraction from the inner region, can redistribute the dark 
matter through a central revirialisation. However the 
most careful study of this process shows the effect to be 
modest, with a reduction in the central halo density by 
at most a factor of two to sb*^. More realistic SPH 
simulations in a cosmological context show that super- 
novae driven turbulent gas motions can impart sufficient 
energy to the dark matter to create a core as large as 400 
parsecs in a Fornax sized galaxy^ 3 -!. This effect requires 
both a significant early central baryon fraction and for 
the Jean's mass to be accurately followed since bulk mo- 
tions are driven by starbursts in giant molecular clouds. 
It will be interesting to compare these experiments with 
higher resolution adaptive mesh techniques including the 
effects of reionisation. 

Over half of disk galaxies have stellar bars which 
can transfer angular momentum to dark mater parti- 
cles through orbital resonances and dynamical friction. 
The magnitude of this process has been debated in the 
literatur o 138 ' 139 ' 140 . However even when a rigid perturber 
mimicking a bar was placed at the center of a CDM 
halo, it only affected the dark matter particles within 
~ 0.001r v ; r ~ 300 pc in our Galaxy. The most recent 
highest resolution study of this process demonstrates that 
the effects of bars on the central dark matter distribution 
is negligible^. 



3.4. Halo shape 

The shapes of dark matter halos can be modified by 
galaxy formatio n 142 ' 143 : Box orbits supporting the tri- 
axial configurations become deformed by the growing 
disbi^ 4 -. Since particles move on eccentric orbits with 
a typical apocentric to pericentric distance of 6:1, halos 
can be visibly affected out to half the virial radius, and 
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become almost spherical close to the galaxy. The change 
in shape depends on the central baryonic fraction, which 
is highest for elliptical galaxies. It is quite low for galaxy 
clusters and dwarf galaxies, their halos should therefore 
barely be affected. The detailed modification of parti- 
cle orbits within the disk region has yet to be explored 
but could also change the predictions for direct detection 
experiments 1 ^. 



3.5. Substructure 

Galaxy formation also leads to the accretion of gas, 
stars and dark matter from satellites into the disk: sys- 
tems on roughly co-planar orbits suffer dynamical friction 
against the disk, which brings them into the disk plane 
where they are disrupte d 145 ' 146 . This process produces 
a dark matter disk, which could contribute a significant 
fraction of the local dark matter density^ 6 -. Satellites on 
more polar orbits on the other hand will pass through the 
disk and experience a quick increase of tidal forces during 
the passage, a so called tidal shock. Such shocks can lead 
to substantial mass loss, their impact depends strongly 
on the orbit of a subhalo and on its mass profil o 74 ' 104 . 

Small, nearby substructures will also lose significant 
amounts of mass as they suffer heating from individual 
star o 147 ' 148 . For the smallest substructures with sizes 
smaller than a few hundred parsecs, impulsive collisional 
heating due to encounters with disk stars dominates their 
mass loss. Over a Hubble time most of their particles will 
be lost into small streams, although often an inner dense 
and strongly bound core is able to survive 1 ^ 8 -. 

The importance of many of these processes remain to 
be quantified which will be an area of activity whilst 
simulators attempt to create realistic galaxies from cos- 
mological initial conditions. They will also play an im- 
portant role in precision cosmology in which future ob- 
servational missions plan to measure the evolution of 
the power spectrum, or the mass distribution on large 
scales to percent level precision, which requires a detailed 
knowledge of how baryons affect the global properties of 
their halos. 



4. PREDICTIONS FOR DIRECT & INDIRECT 
DETECTION EXPERIMENTS 

In this section we will briefly summarize the main pre- 
dictions from cosmological N-body simulations for direct 
and indirect dark matter detection experiments. Some 
of the popular dark matter candidates (e.g. some of the 
WIMPs) would be directly detectable thanks to their tiny 
but just significant enough probability to interact with 
ordinary matter (e.g. atomic nuclei in an underground 
detector). Indirect detection experiments look for 7-rays, 
charged particles or neutrinos produced as two WIMPs 
annihilate each other—. 



4.1. 7-rays from dark matter annihilation 

For the NFW density profile and its variants discussed 
above, the total luminosity from dark matter annihilation 
in a halo scales like 

L cx J p 2 dV 

(7) 

The half light radius is ri/2 — 0.07 ?y mas for the best 
fit profiles in Figure [4} and 0.09 rv m ax for the shallower 
NFW profile. Combining Eqn. [7] with the steep velocity 
functions found for field halos^ 8 - and subhalo s 25 ' 27 ' 31 ' 88 
N(> Vmax) cx dN/dlogV ma x oc V~^ x , implies that small 
field halos and subhalos emit more 7-rays per decade in 
halo size (Vmax or mass) than larger ones. In the lat- 
est galaxy halo simulations all resolved subhalos together 
are about as luminous as the main halo 27 ' 29 i 31 , i.e. the 
boost factor B = Ltotai/^mainhaio is at least two for CDM 
galaxy halos. Extrapolations down to the smallest CDM 
subhalos increase the boost to B = 4 to 16 J—], where 
we assume small mass cut-offs between 10~ 12 and I.OMq 
and we take the nearly constant halo concentrations on 
the smallest scales into account (see 12.2. ip . These boost 
factors imply that small scale structure (< 10 6 M Q ) dom- 
inates the diffuse extragalactic dark matter annihilation 
signal (and also the diffuse galactic signal, except to- 
wards the Galactic cente r 149 ' 150 ). The flux from small 
scale clumps is proportional to the dark matter density 
averaged over larger scales (Figure [6]). 

The spatial and spectral distribution of diffuse gamma 
rays was measured by EGRET— and it is consistent 
with a superposition of (poorly constrained) astrophysi- 
cal diffuse components 1 ^ 2 -. Even under the unrealistic as- 
sumption of a perfect subtraction of astrophysical diffuse 
foregrounds, the DM detection window of the recently 
launched Fermi satellite (formerly know as GLAST) in 
the diffuse galactic component is quite smalli 5 - 3 -. How- 
ever, the significantly improved spatial resolution of 
Fermi relative to EGRET might allow it to detect gamma 
rays from dark matter annihilation in subhalos, both in 
dwarf galaxies and in smaller, dark subhalo s 149 i 150 ' 153 . 
Small scale dumpiness within subhalos increases their 
signal slightly, which makes a somewhat larger number 
of subhalos detectable^. 

4.2. Nearby dark matter distribution and charged 
particles from dark matter annihilation 

Besides 7-rays, dark matter annihilation would pro- 
duce charged particles and anti-particles that, due to 
to magnetic field entanglement, propagate over much 
smaller distances within the Galaxy. Using the local sub- 
halo abundance from Figure [5] and extrapolating down 
to micro-subhalo scales one finds that nearby subhalos 
produce a total flux of 40% of the local smooth halo 
signal—. In other words the local boost factor is 1.4, 
the uncertainty from the extrapolation is about ±0.2. 
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Explaining the positron excess measured by HEAT 154 
and PAMELA 1 — with local dark matter annihilation re- 
quires significantly larger (~ 10 4 ) enhancement o 156 ' 157 . 
When a relatively large subhalo happens to lie within 
1 kpc, the local boost factor increases, but much larger 
values are unlikely: Only 5.2 percent of all random re- 
alizations have a boost factor of 3 or larger (caused by 
a Knax > 3.4kms _1 clump within 1 kpc). In only 1.0 
percent of the cases the boost factor reaches 10 or higher 
due to a nearby, large Kn ax > 5.6 km s -1 subhalo^ 7 -. 



clumps and streams. In pure CDM simulations the rate 
is close to the rate obtained from the standard halo, just 
slightly higher due to the small excess of low velocity 
particles relative to a Gaussian 80 . As for direct detection, 
the presence of accreted dark matter in the Milky Way's 
dis k 146 i 162 is likely the most relevant deviation from the 
standard halo model and further studies are required to 
better understand the properties of dark disks. 



5. SUMMARY AND FUTURE PROSPECTS 



4.3. Local dark matter distribution, direct 
detection and capture in the sun 

Most of the local dark matter is in a smooth 
componen t 27,29 , the probability that the solar system 
currently passes through a subhalo are quite small, 
even when the smallest micro-subhalos are taken into 
accoun t 56 ' 158 . The large number of overlapping streams 
in the inner halo also leads to rather smooth local velocity 
distribution function s 81 ' 95 ' 159 . Even the most prominent 
steams apparent in Figure [7] account for less then one 
percent of the local dark matter density, i.e. even if we 
happened to be located within such a stream today, the 
bulk of detected dark matter particles would still come 
from the "hot" background. Whilst current cosmological 
simulations are able to probe the local density and veloc- 
ity distributions on kpc scales—, alternative methods are 
required to study finer structure s 160 ' 161 . Further stud- 
ies are needed to quantify or exclude the relevance for 
dark matter detection experiments of possible very fine 
grained features in the local six dimensional dark matter 
distribution. 

At 8 kpc from the center of galaxy scale pure CDM 
halos the velocity distributions are peaked (positive kur- 
tosis) , because of the shallower than isothermal potential 
of pure CDM halos fsee 12.2.3)1 . The resulting excess of 
slow and fast particles relative to the Gaussian standard 
halo model, however, is too small to change the interpre- 
tation of direct detection results significantly 80 . Galaxy 
formation likely changed the local velocity distribution 
significantly (c.f. section [3]): For example, the constant 
Milky Way rotation curve implies an isothermal poten- 
tial, which suggest a more Gaussian shape for the real 
velocity distribution of local dark matter particles. 

The shape of the local velocity ellipsoid correlates with 
the shape of the halo: it is radially anisotropic on the 
major axis and tangentially on the minor axis 81 . How- 
ever, while pure CDM halos are elongate d 72 ' 74 , the shape 
of the local Milky Way halo is expected to be fairly 
round 142 ' 143 ' 144 . The dark matter disk (c.f. Section 0§ 
is probably the most drastic deviation from the standard 
halo model, and it may well have significant implications 
for dark matter detectio n 146 ' 162 . 

The time averaged dark matter capture rate relevant 
for neutrino production in the center of the Sun (and 
Earth) 13 smears out the (small) variations due to local 



Numerical simulations provide a reliable method of cal- 
culating the evolution of the dark matter distribution and 
they make robust predictions for the clustering of the 
dark matter component. More recently, but less reliably, 
simulations have begun to resolve the combined growth 
of dark matter and baryonic fluctuations. However there 
is still a lot to learn; the detailed interaction between 
the dark matter and baryons is poorly understood and 
many comparisons with observations and predictions for 
experimental searches rely on extrapolations well below 
the current resolution scale. 

Whilst discrepancies between observations and simula- 
tions have been the subject of much debate in the liter- 
ature, the details of the galaxy formation process needs 
to be resolved in order to fully test the CDM paradigm. 
Whatever the true nature of the dark matter particle 
is, it must not be too different from a cold neutralino 
like particle to maintain all the successes of the model in 
matching large scale structure data and the global prop- 
erties of halos which are mostly in good agreement with 
observations. 

Even as simulations have gone past the "billion par- 
ticle halo" goal, there are still several reasons why it is 
interesting to further increase this resolution. For exam- 
ple, to-date the best resolved CDM satellite capable of 
hosting a dwarf galaxy contains just a few million par- 
ticles and its mass distribution is unresolved on scales 
below a few hundred parsecs. This is the scale within 
which the observations will constrain the dark matter dis- 
tribution through proper motions. Making predictions 
for the phase space structure and substructure within 
the dark matter distribution at the position of the sun 
within the Galaxy is important for direct and indirect 
detection experiments. Resolving the density profile of a 
gravitational collapse as the radius goes to zero is a fas- 
cinating dynamical problem that will guide theoretical 
models towards understanding what processes shape the 
inner structure of dark matter halos. 

Ongoing surveys such as SEGUE 1 ^ 3 . and RAVEiH and 
future surveys such as GAIA and SIM-Lite will provide 
us with detailed measurements of the dynamics and as- 
sembly history of the Galaxy, its stellar halo and satel- 
lite galaxies. The 3-d kinematics of stars from precision 
proper motions will constrain the shape, orientation, den- 
sity law and lumpiness of the Galactic halo. They will 
measure the orbits of the Galactic satellites and the an- 
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gular momentum and orbital anisotropy of the halo stars 
and globular clusters to the outer reaches of the halo. 
The motions of stars in the cores of nearby dSph satel- 
lite galaxies will constrain their phase space density and 
dark matter distributions. Realistic simulations of the 
formation of the Galaxy and its stellar components are 
urgently needed to guide these missions and to aid in 
their interpretation. 

We often hear that this is the precision era of cos- 
mology. Indeed, future precise observations of the mass 
function of dark matter halos^, high redshift 21cm 
tomography^, galaxy surveys such as Pan-STARRS, 
LSST or the Dark Energy Survey and many more, all 
aim to determine the cosmological parameters to a level 
that probes fundamental physics. However there are sev- 
eral problems which arise when numerical simulations are 
used to infer weakly or highly non-linear processes: Dif- 
ferent simulation codes and initial conditions software 
lead to differences in results at the ten percent levels. 
Improving the predictability of simulations by a factor of 
ten will require a significant investment of collaborative 
work. We also need to supply large numbers of large- 
volume cosmological simulations with a high dynamic 
resolution that have numerical and resolution effects fully 
under control. Finally, we are left with large uncertain- 
ties in the role that baryons play in changing some of the 
basic properties of dark matter halos, effects that need 
to be understood before we can proceed in using simu- 



lations to probe fundamental physics rather than testing 
our understanding and modeling of gas-dynamical pro- 
cesses. 

In the past decades and for the foreseeable future, sim- 
ulations of galaxy formation do require guidance from ob- 
servations. In the coming years we will have ground and 
space based telescopes like ALMA or the JWST which 
will reveal details of the galaxy assembly process occur- 
ring at high redshifts. Indeed, one of the ultimate goals 
of simulations is to recreate a single galaxy, or more op- 
timistically, an entire cluster or large volume of well re- 
solved galaxies. The observations are already far ahead 
of the predictive power of numerical simulations, which 
have yet to be able to model the formation of a single pure 
disk galaxy within a cold dark matter halo. Most likely 
this is a problem with resolution, algorithms and sub- 
grid physical processes, which are always present in such 
simulations. More speculatively, this failure could reflect 
a problem with our standard cosmological paradigm. 
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